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ABSTRACT 

With the Navy Prototype Optical Interferometer (NPOI), the binary system 
9 l Orionis C, the most massive member of the Trapezium, was spatially resolved 
over a time period extending from February 2006 to March 2007. The data show 
significant orbital motion over the 14 months, and, after combining the NPOI 
data with previous measurements of the system from the literature, the obser- 
vations span 10 years of the orbit. Our results indicate that the secondary did 
not experience an unusually close periastron passage this year, in contradiction 
to the prediction of a recently published, highly eccentric ~11 year orbit. Fu- 
ture observations of this source will be required to improve the orbital solution. 
Possible implications of the results in terms of system distance are discussed, 
although a main conclusion of this work is that a definitive orbit solution will 
require more time to obtain sufficient phase coverage, and that the interaction 
effects expected at periastron did not occur in 2007. 

Subject headings: binaries: close; open clusters and associations: individual 
(Trapezium); stars: individual {6 1 Orionis C); techniques: interferometric 
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Introduction 



Orion is the nearest example of a giant molecular cloud and is the site of both high- 
mass star formation and a prodigious number of recently-formed stars; the central 2.5 pc 
region of the Orion Neb ula Cluster (ONC ) alone contains ~3500 stars with a combined 



mass exceeding 900 M & ([Hillenbrand 



19971 ). Although star-forming regions such as Taurus 



are closer, these regions are dark clouds ass ociated with low- mass star formation and far 
fewer total stars. Initial star-count studies (ILada et al.l Il99ll ) implied that up to 95% of 
stars formed in clustered environments like Orion, while more recent Spitzer Space Telescop e 
results indicate that at least half of all sta rs originate in dense regions (IMegeath et al.ll2004j ). 
An analysis of binary star distributions (IPatience et al.l |2002| ) suggests that approximately 
70% of field stars may have formed in a clustered environment. Thus, the study of Orion, 
the closest giant star forming cluster, is central to investigating the early history of the 
majority of stars. Furthermore, the variety of cluster morphologies investigated with Spitzer 
observations suggest that OB stars play a signi ficant role in the formation and evolution of 
star formation in clusters (IMegeath et al.l 12004 ) . 



The distance to Orion is a critical parameter that influences the interpretation of many 
of the properties of the region and its members. Measurements range from ~390 pc to ~480 
pc using a variety of methods and assumptions inclu ding observations o f water masers, radio 
sources, an eclipsing b i nary, and statistical analysis dGenzel et al.lll98ll ; IStassun et al.l 12004 



Sandstrom et al.l 120071 ; iMenten et al.l 120071 : Ijeffriesl 120071 . and references therein) . A closer 
distance would imply that the stars are less luminous, and members that are still contracting 
onto the main seque nce are consequently older based on co mparison with theoretical evolu- 



tionary tracks (e.g., ISiess et al.l l2000i ; iPalla fc Stahlerl Il999l ) . Older ages for the stars above 



the main sequence suggest a spread in ages. 

Binary stars yield a model-independent distance with the combination o f a spatially 
resolved orbit and a double-lined radial velocity orbit (e.g., iTorres et al.lll997l ). Given the 
distance to Orion, very high angular resolution is required to separate binary systems which 
exhibit significant orbital motion. The binary system 8 l Ori C (HD 37022) in the Trapezium 
region of Orion has such a close separation tha t it was not detected until speckle observations 
barely resolved the pair (jWeigelt et al.lll999l ) with a separation of less than the diffraction 
limit of the 6m telescope employed in the discovery. The separation has decreased over 
time and is now best monitored with interferometry; a recent orbit fit suggeste d the system 



migh t be just past periastron and completing an orbital cycle within a year (IKraus et al. 



20071 ). In this letter, we present the results of interferometric observations with NPOI that 
add significantly to the binary orbit phase coverage and suggest the orbital period may be 
substantially longer than predicted. 
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The distance to 9 1 Ori C in particular and its phy sical parameters such as mass and 



age are of great importance since the 07 primary (e.g., iGarmany et al.lll982l ). as the most 
massive member of the ONC, has the dominant role in shaping the properties of the sur- 
rounding nebula, and strongly impacts the circumstellar material around the ONC stars. 
'he photoionizing radiation from 8 1 Ori C produces the brightening of many proplyds 

( 



O'Dell et al.lll993l). but also causes the mater ial to escape. Observations of mass loss rates 



(|Johnstone et al 



1998; 



Stoerzer fc Hollenbach 



He nney fc 0'Delllll999l ) and theoretical modeling of the results (e.g., 



19991 ) imply that these stru ctures cannot sur vive for > 10 5 yrs, sub- 



stantially less than the <1— 2 Myr age of the ONC (jHillenbrandlll997l ) . Possible explanations 
for the ap parent contradiction in disk lifetimes and stellar ages include radial orbits for the 
proplyds (jStoerzer fe Hollenbach! Il999l) or a v ery young age for 6 1 Ori C, such as has been 
proposed for 6 1 Ori B (IPalla &: Stahlerl 120011 ). In contrast, recent models of the evolution 
of disk sizes in the ONC ( IClarke 1 120071 ) suggest that disk survival times of 1-2 Myr in the 
uv field produced by 6 l Ori C are possible - consistent with the age of the stellar popula- 
tion. Key to estimating the age of 9 l Ori C is placing the secondary accurately on the H-R 
diagram with a well-measured distance, luminosity, and mass given that the primary has 
already contracted onto the main sequence. We present new magnitude differences which 
will aid in the assessment of the luminosity, but concentrate on the orbital motion which is 
required to estimate the distance. 



2. Observations and Data Reduction 



During Feb. 2006 to Mar. 2007, 9 1 Ori C was observed with NPOI on 6 nights. 
The NPOI array at Anderson Mesa near Flagstaff, AZ, consists of six 50 cm (12 cm clear 
aperture) siderostats w hich can be deployed along a Y-shaped array of vacuum light pipes 
([Armstrong et al.lll998l ). The wavelength coverage spans 550—850 nm over 16 spectral chan- 
nels. Details re gardi ng th e NPOI observation al setup and data recording can be found in 



Hummel et al.l (120031 ) and iBenson et al.l (120031 ) . More recent upgrades relevant to this pro- 
gram are improvements in the acquisition camera sensitivity and calibration of the bias level 
for low count rates; implementation of these two changes enabled the observations of 6 1 Ori 
C. Table 1 provides a log of the NPOI observations. 

The observations of 9 l Ori C were interleaved with e Ori, one of the bright (V=1.70) 
belt stars. The visibilities from e Ori serve to calibrate both variable atmospheric conditions 
and the system response to a point source. The sma ll angular diameter of e Ori - 0.86 ± 
0.16 milli-arcseconds (mas) (IMozurkewich et al.lll99ll ) - and its proximity to 9 l Ori C - less 
than 5° separation - allowed for accurate atmospheric correction. For both the target and 
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calibrator, the fringes were recorded in 2 ms frames for a total scan duration of 30 s before 
switching to the other source. The individual 2 ms data were averaged over a 1 s time period, 
and these 1 s data points were checked for outliers before averaging to generate 30 s averaged 
V 2 measurements. Calibration factors were determined by comparing the observed data of 
e Ori to that expected from a 0.86 mas diameter single star. These calibration factors were 
then applied to observations of 8 1 Ori C. The flagging, averaging, and calibration steps were 
performed with the OYSTER package, as described in Hummel et al. (1998, 2003), except 
that the bias corrections were determined for each star individually. 



3. Results and Analysis 



The calibrated visibilities were fit with a model comprised of two stars with slightly 
resolved stellar surfaces. The primary star diameter was estim ated to be 0.3 mas usin g 
the expected diameter of an 7 star at a distanc e of 450 pc ([Drilling fc Landoltl 120001 ) . 
Assuming the spectral type from lKraus et al.l (120071 ). the secondary diameter was set to half 
that of the primary, rounded to 0.2 mas. The observational setup did not allow us to fit 
for such small diameters convincingly, so we held these values constant. Figure 1 shows 
examples of ca librated V 2 values and the best fit model. Predicted visibilities from a recent 
orbit solution (IKraus et al.l 120071 ) are also plotted. Table 1 lists the epoch, siderostats used, 



number of baselines, maximum baseline length, number of scans, estimated separation and 
position angle with uncertainties, and the position error ellipses. Because each scan yields 
up to 16 V 2 values per baseline and up to 405 V 2 points were obtained during a night, 
individual measurements are not listed, but examples are plotted in Figure 1. Model fits 
to individual baselines include a 180° ambiguity in the position angle, and the values listed 
in Table 1 are chosen based on previous measurements and the inability to fit orbits if the 
companion was located in the opposite quadrant in the 2007 data. Some nights listed in 
Table 1 included closure phase observations, but we defer discussion of these results, which 
have the ability to directly resolve the 180° ambiguity, for a future paper including results 
from calibration binaries. 

While the bandpass does not exactly match the V or R filters, o ur magnitude differ - 
ence measured at 550 nm closely approximates that of the V band ( jZavala et al.l 120071 ). 
Considering all the data, the best estimates for the observed magnitude differences are A 
mag 5 50nm = 1-3 ± 0.3 and A mag7 o n m = 1-2 ± 0.2. The NPOI measurement agrees with 
previous AV estimates from speckle observations t aken at the 6.0 m Special Astrophysical 
Observatory at Mt. Pastukhov in Russia (1.1 mag. IKraus et al.l 120071 ) . 



Orbits based on the previous measurements and the new NPOI data are shown in Figure 
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2. The earliest NPOI data from Fe b. 2006 show only minimal orbital motion from the IOTA 
Dec. 2005 data faraus et al. 2007 ). indicating there is not an offset between the two systems. 
The 2007 NPOI data show significant evolution in the orbit. The measured separations of 
the secondary relative to the primary ar e larger than expe cted from the predicted orbit and 
lag behind the solution in orbital phase (IKraus et al.ll2007l ). The NPOI results indicate that 
a very close periast ron passage d id no t take place this year as suggested by the preliminary 
orbital solutions of IKraus et al.l (120071 ). 



As indicated in Figure 2, only part of the orbit - probably less than one half - is 
covered by the extant data, making any assessment of the orbit fit preliminary. Combining 
all previous position measurements and associated error bars with the NPOI visibilities and 
uncertainties, we fit an orbit using ORBGRID (IHartkopf et al.l Il989l ; iMason et al.l 119991 ) 
and used its solution as a starting point for a least squares solution; uncertainties from 
the covariance matrix of the least squares solution are quoted (Table 2). Both ORBGRID 
and the least squares solution weight the astrometric points in a relative sense, and we set 
these weight s acco rding to the areas of the error ellipses from our Table 1 and Table 3 of 
Kraus et all (hoOTT ). 



We stress that the orbital elements of the best fit based on the current data (Table 2) 
may be modified substantially as more data become available. Compared with the previous 
solution, we find a longer period and a much less eccentric orbit. Our l ower eccentrici ty of 
0.16 is well within the bulk of the distribution for T Tauri binaries (e.g.. lMathieulll9941 ). the 



lower mass counterparts of 6 1 Ori C. In contrast, the earlier solution found an extremely 
high eccentricity of 0.91—0.93 (Kraus et al. 2007). 

Although it is premature to calculate a robust distance to the Orion Nebula Cluster 
from our data, if we assume a total system mass of 40 M the orbital elements in Table [2] 
give a dynamical parallax distance of 730 pc — unr ealistically large con sidering the distance 
to the background high extinction molecular gas (IGenzel et al.l 119811). The 40 M va lue 



is the minimum mass estimate from th e evolutionary track s of IWalborn fc Nichols! (119941 ). 
using a T eff for an 07 star of 36,000 K (IMassey et al.ll2005l ). The values of Amag 5 5o n m and 
Amag 70 onm b etween the primary and secondary stars imply at the latest a B2 secondary 
spectral type ([Drilling fc Landoltl 120001 ). with a corresponding T e g of 28,000 K. Thus, 40 M 
is a lower limit; larger masses would imply an even greater distance. 

Figure 3 summarizes our dynamical parallax measurement. The uncertainty in the value 
of a 3 /P 2 (Table [2]), directly related to parallax, does not yield a reliable estimate of distance 
at this time. Increasing the semi-major axis by 1 a and reducing the period by 1 a yields 
a low value for the distance of 344 pc for the same total mass. Clearly, the available data 
do not significantly constrain the distance to 9 1 Ori C. Figure 3 compares the best-fit semi- 
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major axis and period, with their associated 1 a and 2 a error ellipses, to the values expected 
for the two Orion distances estimated from the radio star and the maser. To explore fits to 
our data with more physical distances imposed, we selected a range of nine periods from 10 
to 26 years and determined the corresponding semi-major axes for distances of 390 pc and 
480 pc. For both distances, fits were obtained for periods of ~22 years with orbital elements 
that agreed to within ~ 3a of the best fit elements (Table [2]). Fits for other periods were 
significantly (by many sigma) worse. 

Further monitoring of the orbit of 6 1 Ori C is required to decrease the errors in the 
orbital elements and provide a reliable dynamical distance. Given the proximity to periastron 
passage, continued observations are particularly important over the next months and years. 
Resolving the distance and mass of 8 1 Ori C, and revealing the nature of its interactions 
with the local environment, will provide important insight into the closest region of massive 
star formation. 

The Navy Prototype Optical Interferometer is a joint project of the Naval Research 
Laboratory and the US Naval Observatory, in cooperation with Lowell Observatory, and is 
funded by the Office of Naval Research and the Oceanographer of the Navy. The authors 
would like to thank Jim Benson and the NPOI observational support staff Dave Allen, Jim 
Clark, Brit O'Neill, Susan Strosahl, Dale Theiling, Josh Walton and Ron Winner whose 
efforts made this project possible. We thank Phil Massey, Nathan Mayne for helpful discus- 
sions, and Nat White for assistance with observing arrangements. RTZ thanks the Michelson 
Science Center for an invitation to visit which helped initiate this work, and JP gratefully 
acknowledges funding from the Michelson Fellowship Program. 
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Table 1. NPOI Observations and V 2 Model Fit Results 



UT Date 


Julian Year 


Siderostats 


# b.l. 


Max b.l. 




P 




e 








"4, (°) 


(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


(7) 


(8) 


(9) 


(10) 


(ii) 


(12) 


(13) 


2006 Feb 24 


2006.1486 


AC-AE-AW 


2 


22.2 


171 


11.80 


1.11 


152.3 


3.5 


1.20 


0.54 


178.0 


2007 Feb 22 


2007.1425 


AE-AW-W7 


2 


37.4 


118 


11.94 


0.31 


88.1 


5.2 


1.09 


0.28 


170.4 


2007 Feb 25 


2007.1507 


AE-AW-W7 


2 


37.4 


60 


12.13 


1.58 


92.9 


8.8 


2.41 


0.39 


142.8 


2007 Mar 06 


2007.1753 


AE-AW-AN-W7 


4 


38.1 


405 


12.17 


0.37 


86.6 


2.1 


0.46 


0.36 


157.6 


2007 Mar 17 


2007.2055 


AE-AW-AN-W7 


4 


38.1 


135 


12.28 


0.41 


85.6 


1.9 


0.46 


0.35 


42.2 


2007 Mar 20 


2007.2137 


AE-AW-AN-W7 


4 


38.1 


200 


12.14 


0.43 


83.0 


2.3 


0.49 


0.42 


158.2 



Note. — Col. (1): UT date of NPOI observation Col. (2): Julian Year of NPOI observation. Col. (3): Siderostats used Col. (4): Number 
of independent baselines Col. (5): Max. baseline length (m) Col. (6): Number of V 2 measurements Col. (7): Fitted binary separation (mas) 
Col. (8): Error in p (mas) Col. (9): Fitted binary position angle (°).; this angle is chosen as it is a smooth extension of previous results and no 
orbital solutions could be found using 6 + 180° Col. (10): Error in 6 (°) Col. (11): Semi-major axis of error ellipse (mas). Col. (12): Semi-minor 
axis of error ellipse (mas) Col. 13: Position angle of error ellipse. 
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Table 2. Orbital Elements 



Data 


Value 


a (mas) 


41 ± 14 


i (deg) 


107.2 ± 3.5 


n (deg) 


208.8 ± 3.7 


e 


0.16 ± 0.14 


u (deg) 


96.9 ± 118.7 


T (JY) 


2007.0 ± 5.9 


T (JD) 


2454101 ± 2154 


P (days) 


9497 ± 1461 


P (years) 


26 ± 13 


a 3 /P 2 (mas 3 /yr 2 ) 


103 ± 146 


X 2 /dof 


2.6xl0- 6 
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Fig. 1. — Left: Panels showing 3 of 9 scans of the calibrated V 2 of 6 1 Ori C observed with the 
AE-AW baseline on 2007 Mar 06 overlaid with a model of separation 12.17 mas and position 
angle 86.6° (Table 1). Right: A set of three panels shows the same data but o verlaid using 
the pr edicted separation of 8.6 mas at a position angle of 86.28° for orbit 1 of iKraus et al. 
(120071 ) on that date. The errors include both a statistical error and an estimate of the 
uncertainty of the calibration. 
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Fig. 2. — The best-fit orbit of 6 1 Ori C with orbital elements listed in Table [2] is plotted as 
a solid line. This fit inc orporates th e NPO I astrometric results from Table 1 and previous 
measurements given in iKraus et al.l (120071 ). Error ellipses for the astrometric points are 
shown along with a vector i ndicating the peri astron point. The dashed line shows the 
predicted orbit (Orbitl) from lKraus et al.l (120071 ). We display the closely spaced 2007 NPOI 
astrometric solutions in the inset for clarity. The next set of observations should demonstrate 
decisively the validity of a longer period solution. 




20 40 

Period (years) 

Fig. 3. — The best-fit period and semi-major axis from Table 2 are plotted with a cross, 
and the superimposed error ellipses (solid lines) represent the 1 and 2 a errors of the orbital 
elements in Table 2. For a total system mass of 40 M Q , the period and semi-major axis 
values consistent with distances of 480 pc (dashed) and 380 pc (dotted) are also indicated. 
At the ~1 a level, the current orbit fit is consistent with essentially all measured distances 
to 9 1 Ori C. Given the limited phase coverage of the orbit a definitive distance measurement 
will require continued observations. 



